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Aqueous humor (AH) is one of the body fluids in the eye, which is known to be related with various ocular 
diseases, but the complete RNAs characteristic of the AH in patients is not yet known. The aim of this study 
was, with a microarray analysis, to reveal the disease-related extracellular miRNAs profiles in individual 
patients AH. 100 ul of AH was collected by anterior chamber paracentesis from 10 glaucoma, 5 cataract, and 
5 epiretinal membrane patients. The extracted total RNAs were shorter than 200 nt, and their amount was 
5.27 ± 0.41 ng in average. Among 530.5 ± 44.6 miRNA types detected in each sample with a microarray 
detectable 2019 types of matured miRNAs, 172 miRNAs were detected in all 10 glaucoma or control 
patients. From the glaucoma group, 11 significantly up-regulated and 18 significantly down-regulated 
miRNAs (P < 0.05 for both) were found to have areas under the curve better than 0.74 in a receiver operating 
characteristic analysis. They also formed a cluster composed only of glaucoma patients in a hierarchal cluster 
analysis. AH had a possibility of becoming a source of miRNA that can serve as a biomarker and a 
therapeutic target. 



I ncreasing the quality of care expected by patients requires the development of improved diagnostic, prognostic 
I and therapeutic monitoring tools for disease. Biomarkers of disease, including extracellular microRNAs 

I (miRNAs) in bodily fluids, may play a role in achieving this purpose, as well as elucidating pathogenic 
mechanisms. MiRNAs have recently attracted attention for the prospect of earlier diagnosis and more precise 
application of existing therapies 1 . 

MiRNAs are small, non-coding RNA molecules, consisting of 19-22 nucleotides. They have been intensely 
studied since the discovery two decades ago of the role of the lin-4 gene in regulating protein abundance 2,3 with the 
general goals of understanding their role in translational repression and messenger RNA cleavage in various 
biological processes 4 6 , and elucidating pathologies in these processes 6-9 . After the recent discovery of circulating 
extracellular miRNAs in bodily fluids (e.g. plasma and serum), new searches for disease-related extracellular 
miRNAs have started in various fields 4,10-13 . 

There are three sources of extracellular miRNAs in the eye: tears, the vitreous humor, and the aqueous humor 
(AH) 11,14,15 . Studies of ocular fluids and other cellular tissues that used a proteomic approach have identified 
several disease-related proteins 1619 including multiple Alzheimer's disease related peptides in the AH of eyes with 
glaucoma 20,21 . In comparison, research into ocular extracellular miRNAs is only at an early stage 11,14,15 and still has 
great potential to uncover novel, accurate and valuable biomarkers and therapeutic targets. 

Glaucoma is one of the leading causes of blindness worldwide 22 24 . The axons of retinal ganglion cells, which 
extend to the brain through the optic nerve, are irreversibly damaged in glaucoma, but with early detection and 
proper medical treatment, it is possible to slow the progress of this damage. The present study focused on finding 
AH miRNAs related to glaucoma. The AH is the most frequently sampled ocular fluid in examinations and 
clinical studies of ocular disease. It is secreted by the ciliary epithelium, which supplies nutrients and removes 
metabolic waste from the avascular tissues of the eye such as the lens and cornea. The AH is excreted through the 
angle of the anterior chamber, which consists of the trabecular meshwork, uvea, and sclera. Pressure regulation 
and dysregulation in this system is associated with glaucoma 25 . Luna et al reported that miR-29b has a potential 
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role in glaucoma with an in vitro assay of the trabecular meshwork 
cells 26,27 . MiRNAs in the AH may have a key role in pathological 
conditions related to glaucoma, but the entire miRNA profile of 
the glaucomatous AH has not yet been determined. 

The discovery of novel miRNAs has been accelerating with 
advances in analytical technologies 28 , and the total number of known 
human mature miRNAs has reached 2019 (miRBase Version 19, 
released on Aug. 1, 2012). The possibility of finding new and more 
useful biomarkers thus continues to rise. In this study, we used a 
microarray system able to detect all 2019 known mature miRNAs to 
analyze individual AH samples from patients with glaucoma, catar- 
acts, and epiretinal membrane (ERM), with the aim of revealing 
specific diseases, individual expression characteristics, and the num- 
ber of detectible miRNA types in the AH. 

Results 

Total RNA extraction from AH of ocular patients. Table 1 and 
supplementary Table SI shows patient characteristics, quantity of 
RNA yielded, and the number of miRNAs detected by the 
microarray system. There was no significant difference in age 
between the 10 glaucoma patients (72.8 ± 3.0 years old) and the 
10 control patients (cataract and epiretinal membrane, 71.1 ± 2.6 
years old), who participated in this study (Supplementary Figure. S2a 
online). Intra ocular pressure of glaucoma group was significantly 
higher than control group (Supplementary Figure. S2b online). We 
were able to obtain high quality purified total RNA with the 
bioanalyzer. There was a significant quantity of RNA shorter than 
200 nt, and no ribosomal RNA longer than 200 nt (specifically, 18S 
and28S rRNA were absent) (Figure la). From a 100 ul sample of AH, 
we obtained 5.27 ± 0.54 ng of purified total RNA from the glaucoma 
patients and 6.19 ± 0.62 ng from the 10 control patients (the average 
for both groups was 5.73 ± 0.41 ng). There was no significant 
difference in the quantity of purified total RNA obtained from the 
two groups (Figure lb). 

Microarray analysis of miRNAs in AH. Even with the small amount 
of purified total RNA and without any amplification steps, bright 



spots were present on the microarray (Figure 2a). An average of 530.5 
± 44.6 miRNA types were detected in the 20 samples, with no 
correlation to the quantity of purified total RNA (R 2 = 0.0197, 
Table 1, Supplementary Fig. S2 online), and no significant 
difference in the number of miRNAs detected in the glaucoma 
patients (500.0 ± 59.7) or control patients (561.0 ± 68.2) (Figure 2b). 

Statistical analysis of microarray data of AH miRNAs. All of the 

values observed were normalized to the expression value of hsa-miR- 
3940-5p, which was detected in all the samples (p = 0.985, ratio of 
1.005, Supplementary Table S2 online). 57 miRNAs showed a 
statistically significant difference (p < 0.05, q < 0.2) in expression 
levels in the control patients and glaucoma patients (Supplementary 
Figure S2 online). In this group, hsa-let-7b-3p (Supplementary Table 
S3 online), has-miR-4507, has-miR-3620-5p, has-miR-1587, and 
has-miR-4484 were most significantly different (p < 0.01). 

(a) (b) 




Figure 1 | (a) Total RNA validation, screens capture of Agilent 2100 
Bioanalyzer electropherograms of representative human AH samples 
(total RNA), (b) comparison of total RNA quantity in glaucoma and 
control patients. 
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Figure 2 | (a) MiRNA microarray analysis, representative scanned image 
of a microarray capable of detecting 2019 miRNAs, and (b) comparison 
of detected miRNAs number in glaucoma and control patients. 



The numbers of miRNA types detected in the all of 10 glaucoma 
patients and the all of 10 control patients were 158 and 135, respect- 
ively (Figure 3a). 121 of these miRNAs (group A) were detected in all 
20 patients. 14 miRNAs (group B) were detected in all glaucoma 
patients, but not detected in any control patients. 37 miRNAs (group 



C) were detected in all control patients, but not detected in any 
glaucoma patients. Among the 121 miRNAs of group A, 18 
miRNAs showed a statistically significant difference (<0.05) in 
expression level in the control patients and glaucoma patients, as 
shown in Fig. 3b and supplementary Table S4 online. Among these 
18 miRNAs, 8 were up-regulated and 10 were down- regulated in the 
glaucoma patients. Among the 14 miRNAs of group B, 3 miRNAs 
were up-regulated in the glaucoma patients, a significant difference 
(<0.05) (Figure 3c, and Supplementary Table S5 online). Among the 
37 miRNAs of group C, 8 miRNAs were down-regulated in the 
glaucoma patients, a significant difference (<0.05) dFigure 3c, and 
Supplementary Table S6 online). A receiver operating characteristic 
(ROC) analysis of the 11 up-regulated and 18 down-regulated 
miRNAs showed that they all had an area under the curve (AUC) 
greater than 0.74 (Figure 3b, c, and d). A hierarchal cluster analysis 
with these 29 miRNA markers showed a cluster composed only of 
glaucoma patients was formed (Figure 4). 

Targets and pathway analysis of glaucoma specific AH miRNAs. 

DIANA-microT analysis was performed on 29 miRNA markers and 
hsa-let-7b-3p. Sixty-two and 266 sets of up- and down- regulated 
miRNAs/targets were predicted to have miTG scores > 0.97000 
(Supplementary Table S7 and S8 online). Ten and 9 molecules, 
respectively, were targeted by more than one up- or down- 
regulated miRNA (Table 2). The top 2 molecular networks 
predicted by the IP A analysis were (1) RNA post-transcriptional 
modification, developmental disorders, hereditary disorders, and 



(a) 



158 miRNAs detected 
in all 10 control patients 



miRNAs detected in all 10 control patients 
not detected in all glaucoma patients 
(Group C) 



(b) Group A 




135 miRNAs detected 
in all 10 glaucoma patients 



miRNAs detected in all 10 glaucoma patients 
not detected in all control patients 

miRNAs detected in all 20 patients (Group B) 

(Group A) 

(c) Group B 



Name of miRNA 


Average 
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p-value Fold 
(t-test) change 


Difference 
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0.095 


0.069 
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* AUC = area under the curve. 



Figure 3 | (a) Venn diagram of the miRNAs detected in all control patients and all glaucoma patients. We detected 121 miRNAs in all 20 patients (group 
A), 14 miRNAs in all glaucoma patients, but not in all control patients (group B), and 37 miRNAs in all control patients, but not in all glaucoma patients 
(group C). (b) List of 18 miRNAs with a significant difference in glaucoma and control patients, screened from group A. (c) List of 3 miRNAs up- 
regulated in glaucoma patients and exhibiting a significant difference in the glaucoma and control patients, screened from group B. (d) List of 8 miRNAs 
down-regulated in glaucoma patients and exhibiting a significant difference in glaucoma and control patients, screened from group C. 
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Figure 4 | Cluster analysis of selected 29 miRNA markers for glaucoma. 



(2) tissue development, neurological diseases, and auditory diseases 
(Table 3). The top 2 predicted canonical pathways were (1) protein 
kinase A signaling and (2) calcium signaling (Table 4). 

Discussion 

The AH is an attractive source of novel miRNA biomarkers of ocular 
diseases, because of its accessibility, specificity, independence from 
other organs and novelty. Even though many researchers have 
reported various risk factors and biomarkers for ocular diseases epi- 
demiologically and genetically 29 35 , studies of AH RNA are rare. The 
primary reason is that the AH is not a cellular sample, and the 
amount of total RNA in the AH is so small that it is undetectable 
with the usual spectophotometrical procedures. Even against such a 



challenging background, Dunmire JJ et al. reported a preliminary 
result suggesting that miRNAs exist in the AH in 2013 14 , which raised 
interest in further study of AH RNA profiles and their use as ocular 
disease biomarkers. 

In the present study, we evaluated the concentration of total RNA 
in the AH (14-85 ng/ml) for the first time, and found that it is similar 
to human breast milk (9.7-228.2 ng/ml) as reported by Kosaka N. 
et al. 10 , but less than in other body fluids (e.g. amniotic fluid, tears, 
seminal fluids, and so on) as reported by Jessica AW 11 . Addi- 
tionally, the size distribution of RNA in the AH (<200 nt, Fig. la) 
is shorter than in breast milk (OOO nt). The high miRN A/small 
RNA purity in the AH may compensate for miRNAs detection and 
relevancy as a source of biomarkers. 
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One of the advancements of our study was the larger number of 
tested miRNAs than in any other study of miRNAs in other body 
fluids (e.g. 264 targets with a PCR array for an AH mixture from 5 
cataract patients 14 , 723 targets with a micro array for breast milk 10 , 
and 706 targets with a PCR array for 12 body fluids 11 ). Detected 
miRNAs in an individual AH sample reached five hundred on aver- 
age, which was approx. 5 times more than that were detected in a 
mixture of 5 AH samples. Furthermore, our study is the first to reveal 
that individual miRNA profiles vary for every patient, and that the 
number of miRNA(s) detected in common was limited to 121 (26%, 
in average) in 20 patients. 

This is the first report to identify a number of miRNAs in the AH, 
and compare them between glaucoma and control patients (cataract, 
and ERM) (Figure 3b, c, and d). We also found a higher AUC for 
these miRNAs in the glaucoma patients. In addition, Hsa-let-7b-3p 
was revealed as having the best odds ratio for glaucoma (7 glaucoma 
patients and 0 control patients) in this study (Supplementary Table 
S2 online). Furthermore, combining these markers has the potential 
to increase the sensitivity of glaucoma diagnosis (Figure 4). 

We also found that miRNA profiles had the potential to reveal the 
stage of pathology, which is not possible with a genome analysis. For 
example, Let-7b miRNA, which had the best odds ratio for glaucoma 
in this study, has been reported to be an age-related marker of catar- 
acts when expressed in the lens epithelium 36 , as well as a down- 
regulation marker of tumorigenesis and retinoblastomas 37 . 
Additionally, bioinformatics predicted which molecular targets 
(Table 2, Supplementary Table S6 online), molecular networks 
(Table 3), and canonical pathways (Table 4) were related to the 
detected miRNA biomarker candidates. Several predicted target 
molecules such as ENFB3, ESRRG, and BCL2 family (Table 2, 
Supplementary Table S7 online) have been already reported as glauc- 
oma-related-molecules in clinical study and/or animal model 38 40 . 



The predicted data obtained in this study, may lead to useful findings 
on the pathology of glaucoma and may reveal novel therapeutic 
targets. 

Glaucoma can be categorized into primary open angle glaucoma 
(POAG), normal tension glaucoma (NTG), primary angle closure 
glaucoma (PACG), and other subtypes, and their pathological con- 
ditions are extremely varied. In addition, medication for patients has 
a potential to influence in AH miRNA profiles. In a future study, we 
expect to find miRNA biomarkers that can diagnose these glaucoma 
subtypes, pathological conditions, and pharmacological effects. We 
also expect that discovery of new miRNAs will be dramatically accel- 
erated by the development of new analyzer machines, such as next- 
generation sequencers 28,4143 . Exhaustive analysis of miRNAs is 
expected to be carried out periodically to update our knowledge of 
biomarkers and therapeutic targets. 

In summery, we found short RNA in the AH of our patients, and 
profiling with a microarray system identified several candidate 
miRNAs related to glaucoma. Additionally, more than a hundred 
miRNAs were found in most individual AH miRNA profiles that had 
potential as biomarkers and therapeutic targets for various ocular 
and other diseases. 

Methods 

AH sampling. This research (University Hospital Medical Information Network; 
UMIN Study ID N.: UMIN000011121) was conducted in accordance with the 
Declaration of Helsinki and Tohoku University Medical School's Institutional Review 
Board approved the protocol (2012-1-546). Informed consent was obtained from 
patients undergoing glaucoma, cataract, and epiretinal membrane surgery at Tohoku 
University Hospital, Sendai, Miyagi, Japan. AH was collected from patients who met 
the inclusion criterion of having had no history of cancer, asthma, and ocular diseases 
other than glaucoma, cataract, and epiretinal membrane. 

The enrolment and grouping criteria were listed as follows: (1) POAG was diag- 
nosed as glaucomatous optic neuropathy on fundoscopic examination, characteristic 
visual field defects on the Humphrey Field Analyzer according to the Anderson- 



Table 3 | Top 5 molecular networks predicted by an IPA analysis of miRNA targets 

Score Focus Molecules Top Diseases and Functions 



48 

37 
33 
32 
30 



27 RNA Posr-Transcriptional Modification, Developmental Disorder, Hereditary Disorder 

22 Tissue Development, Neurological Disease, Auditory Disease 

20 Connective Tissue Disorders, Dermatological Diseases and Conditions, Gastrointestinal Disease 

22 Organ Morphology, Renal and Urological System Development and Function, Cellular Assembly and Organization 

1 9 Cellular Assembly and Organization, DNA Replication, Recombination, and Repair, Cell Morphology 
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Table 4 | Top 5 canonical pathways predicted by an IPA analysis of miRNA targets 
Ingenuity Canonical Pathways "'°9(p 

-value) R a tio 



Protein Kinase A Signaling 1 .60E+00 2.46E-02 

Calcium Signaling ' 1 .54E+00 2.76E-02 

Relaxin Signaling 1 .50E+00 3.05E-02 

Clathrin-mediatecl Endocytosis Signaling 1 .49E+00 3.03E-02 

Role of CHK Proteins in Cell Cycle ' 1 .47E+00 5.08E-02 
Checkpoint Control 



Molecules 



PDE2A,NFAT5,FLNA,DUSP1 ,GNG2,CHP1 ,ELK1 ,PTPN1 2,CDC25A,PDE1 C 

NFAT5,CHRFAM7A,CHP1,CHRNA7JPM4,TRPC3 

PDE2A,GNAO! ,GNG2,ELK1 ,PDE1 C 

SNX9,RAB5A,USP9X,CHP1,NUMB,MY01E 

PPM1J,CDKN1A,CDC25A 



Pattela classification confirmed in at least two visual field examinations, open 
anterior -chamber angles on gonioscopy, and elevated intraocular pressure (intra- 
ocular pressure (IOP) > 21 mm Hg by the Goldmann applanation tonometer); (2) 
The diagnosis criteria for NTG were identical to those for POAG, except that the IOP 
never exceeded 21 mm Hg; (3) PACG was defined as glaucomatous optic neuropathy, 
above characteristic visual field defect, and an occludible angle; (4) Pseudoexfoliation 
glaucoma (PEX) was defined in this study as glaucomatous optic neuropathy with the 
same characteristic visual field defects as the previously mentioned forms of glauc- 
oma, an open anterior chamber angle in a gonioscopic examination, and the finding 
of characteristic exfoliative material on the anterior lens surface and/or in the iris in a 
slit-lamp examination, in one or both eyes. The presence of abnormal visual field 
defects meant that the results of a glaucoma hemifield test were outside normal limits 
and that a cluster of three or more nonedge points were present, all depressed on the 
pattern deviation plot at P < 0.05 (with one depressed at P < 0.01), as well as that the 
corrected pattern SD was significant at P < 0.05. An occludible angle was classified 
when the posterior trabecular meshwork could not be seen over an angle of 180 
degrees or more without indentation. 

Cataract and macular disorders such ERM are the most common non-exudative 
retinal diseases used as control groups in clinical ophthalmological studies, including 
those of glaucoma 44 . The eyes with ERM included in this study were exclusively eyes 
with idiopathic ERM, and did not have any other macular abnormalities. Eyes with 
secondary ERM (e.g., attributable to diabetic retinopathy, venous occlusion, retinal 
detachment, uveitis, or trauma) were thus excluded. 

AH was collected from patients who met the inclusion criterion of having had no 
history of cancer, asthma, and ocular diseases other than glaucoma, cataract, and 
epiretinal membrane. 

AH sampling and RNA isolation. Approximately 100 ul of AH was collected from 
each patient by anterior chamber paracentesis, using a 30-gauge needle inserted 
through the peripheral cornea at the beginning of the procedure. The needle did not 
contact any iris or lens tissue during the sample collection. AH was transferred to 
1.5 ml siliconized tubes and immediately placed on ice. Samples were subsequently 
centrifuged for 10 min at 300 g for removing cellular components, and its aqueous 
phase was gently collected. The AH was mixed with 700 ul QIAzol Lysis reagent and 
stored at — 80°C until further processing. Purification of miRNA from thawed AH/ 
QIAzol samples were performed with miRNeasy Mini Kit (Qiagen, Valencia, CA). 
And fraction of small RNAs (<200 nucleotides) was conformed by Agilent 2100 
bioanalyzer (Agilent teqnology) with a RNA 6000 Pico kit. 

miRNA array preparation and analysis. MiRNA in purified RNA was measured 
using the Toray Industries miRNA analysis system, in which AH miRNA samples 
were hybridized to 3D-Gene human miRNA ver. 1.60 chips containing 2019 miRNAs 
(Toray Industries, Inc., Tokyo, Japan). MiRNA gene expression data were scaled by 
global normalization, and differential expression was analyzed 45 . Two-tailed, Welch's 
t-test was used to detect significant associations. P-values were adjusted for multiple 
testing based on the false discovery rate with the Bioconductor package q-value 45 . 

Bioinformatical analysis. After a determination was made of which miRNAs had 
significantly different levels of expression in the control subjects and glaucoma 
patients, the molecular targets of these miRNAs were predicted with DIANA-microT 
v3.0 (http://diana.cslab.ece.ntua.gr/microT/) 47 . Pathway and global functional 
analyses were performed with IPA software 48,49 . The dataset uploaded to the IPA 
application was the list of target molecules with miTG scores > 0.97000 from the 
DIANA-microT prediction. Values of plus one and minus one were assigned to the 
molecules predicted to be targeted by down- and up- regulated miRNAs, respectively. 
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